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ABSTRACT 

A  single-scattering  atmospheric  radiative  transfer  model  is  de¬ 
rived  and  described.  This  model,  based  only  on  ground  level  measure¬ 
ments,  can  be  used  to  predict  the  spectral  and  directional  character¬ 
istics  of  radiation  incident  cr.  a  ground  target  from  the  sun  and  the 
cloudless  daytime  sky  in  the  absence  of  discontinuities  such  as  inversion 
layers.  Wavelengtns  from  0.4  to  1.7  microns  are  considered. 

The  model  makes  use  of  available  information  on  the  scattering 
and  absorbing  properties  of  atmospheric  constituents  and  the  distri¬ 
bution  of  t.iese  constituents  with  altitude.  The  final  result  is  in 
the  form  of  a  function  that  can  be  integrated  by  computer  to  determine 
the  spectral  irradiance  on  the  target  from  a  preselected  segment  of 
the  clear  daytime  sky.  Detailed  discussions  of  error  sources  and  model 
limitations  are  included. 


3 


Preceding  page  blank 


TABLE  OF  CONTENTS 

Page 

ABSTRACT .  3 

LIST  OF  ILLUSTRATIONS .  7 

I.  BACKGROUND .  9 

II.  PROBLEM  DESCRIPTION  .  10 

III.  EVALUATION  OF  MODEL  INPUT  VARIABLES  .  12 

A.  Evaluation  of  Hg . 12 

B.  Evaluation  of  fldV  .  12 

C.  Evaluation  of  B(5) . 13 

D.  Evaluation  of  T^  . 18 

IV.  UTILIZATION  OF  THE  MODEL . 26 

V.  CONCLUSIONS . 29 

REFERENCES . 32 

APPENDIX  A  -  ILLUSTRATIONS . 35 

DISTRIBUTION  LIST . 61 


5 

Preceding  page  Hank 


ikwZS&TfS' 


LIST  OF  ILLUSTRATIONS 

Figure  Page 

1.  Format  Used  in  Predicting  Sky  Radiance- . 36 

2.  H$  as  a  Function  of  Wavelength . 37 

3.  6  as  a  Function  of  6.  and  <f>  for  6  =  20  Degrees . 38 

L  S 

4.  6  as  a  Function  of  6.  and  <j>  for  0  =  40  Degrees . 39 

L  S 

5.  6  as  a  Function  of  6.  and  $  for  0  =  60  Degrees . 40 

L  S 

6.  6  as  a  Function  of  0,  and  <f>  for  0  =  80  Degrees . 4] 

7.  The  Scattering  Area  Ratio  as  a  Function  of  2ttR/A  ....  42 

8.  Relative  Wove  length  Response  of  the  Human  Eye . 43 

9.  Droplet  Concentration  as  a  Function  of  the  Visibility  .  .  44 


10.  Mie  Scattering  Coefficient  as  a  Function  of  the 

Visibility  and  Wavelength  .  45 

11.  Normalized  Directional  Mie  Scattering  Coefficient  ....  46 

12.  Altitude  Variation  of  the  Visual  Spectrum  Volume 

Scattering  Coefficient  after  Waldram14 . 47 

13.  G(A)  as  a  Function  of  Wavelength . 48 

14.  Mie  Transmission  Geometry  for  V  >  5(A) . 49 

15.  Carbon  Dioxide  Optical  Slant  Range  as  a  Function 

of  8  ,  0.  ,  and  r . 80 

S  Li 

16.  Carbon  Dioxide  Optical  Slant  Range  as  a  Function 

of  0  ,  9.  ,  and  . . 51 

S  L 

17.  Carbon  Dioxide  Optical  Slant  Range  as  a  Function 

of  0_,  0, ,  and  . 52 

5  Li 

18.  Water  Vapor  Concentration  as  a  Function  of  Altitude  ...  S3 

19.  Water  Vapor  Concentration  as  a  Function  of 

Meteorological  Conditions  .  54 


7 


fnaiist  pags  Hank 


Page 


LIST  OF  ILLUSTRATIONS  (CONT.) 


Figure 


20.  Water  Vapor  Optical  Slant  Range  as  a  Function  of 

0  ,  0,  ,  and  r .  55 

S 

21.  Water  Vapor  Optical  Slant  Range  as  a  Function  of 

6,0.,  and  r . 56 

S  L 

22.  Water  Vapor  Optical  Slant  Range  as  a  Function  of 

0,6.,  and  r . 57 

S  L 


23.  Mie  Transmission  Geometry  for  V  <  G(A)  and  r  <  1/sin  0^  .  58 

24.  Mie  Transmission  Geometry  for  V  <  G(X)  and  r  >  1/sin  0^  .  59 

25.  Altitude  Variation  of  the  Visual  Spectrum  Mie 

Volume  Scattering  Coefficient  as  Predicted  by 

the  Model . 60 


8 


I .  BACKGROUND 

Reflectance  modeling  for  natural  sources  of  irradiation  is  a  part 
of  a  comprehensive  ground  target  signature  modeling  program.  The 
objective  of  this  research  is  to  develop  a  unified  methodology,  vis. 
predictive  and  empirical  modeling  of  the  basic  fundamental  physical 
phenomena  that  are  involved  in  the  application  of  ground  target  and 
background  signatures  to  the  design  and  evaluation  of  target  acquisi¬ 
tion  systems,  terminal  homing  systems  and  countermeasures. 

The  natural  sources  of  irradiation  to  be  considered  include  the 
sun,  s\y  and  ground  during  the  day  and  the  moon,  sky  and  ground  at 
night.  The  precedure  to  be  followed  in  modeling  these  sources  is  the 
same  for  the  daytime  and  nighttime  cases  except  that  at  night  an  air- 
glow  term  must  be  added  to  the  sky  input.  The  wavelengths  of  interest 
start  in  the  visual  part  of  the  spectrum  CO .4  micron)  and  extend  to 
the  wavelength  at  which  reflected  energy  is  negligible  compared  to  the 
emitted  energy.  This  upper  wavelength  has  been  selected  to  be  3.5 
microns.  The  wavelength  increment  must  be  small  enough  to  permit  sen¬ 
sor  optimization  studies  to  be  performed  but  large  enough  that  the 
data  handling  procedure  does  not  become  too  complicated.  The  wave¬ 
length  increment  selected  for  this  program  is  0.05  micron.  For 
simplicity,  only  single  scattering  has  been  considered.  The  effect 
of  multiple  scattering  can  be  large  in  some  cases  and  will  therefore 
be  addressed  later  in  the  program. 

Once  the  sources  are  adequately  charactei ized,  the  radiation- 
target  interaction  must  be  considered.  The  only  phenomenon  treated 
in  this  portion  of  the  program  is  reflection.  It  is  necessary'  to 
know  the  bidirectional  reflection  characteristics  of  target  components 
as  a  function  :f  wavelength.  In  order  to  specify  the  angles  of  inci¬ 
dence  and  reflection,  it  is  necessary  to  know  the  directions  of  the 
sources,  the  direction  of  the  sensor,  and  the  direction  of  the  normal 
to  the  target  surface  at  the  point  in  question.  The  direction  of  the 
sensor  is  an  assigned  parameter.  The  target  surface  normal  can  be 
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defined  based  on  a  combinatorial  geometry  algorithm  developed  by 

1* 

Armament  Systems,  Inc.  This  report  describes  an  initial  approach 
to  modeling  the  direction,  intensity,  and  spectral  characteristics  of 
each  of  the  relevant  sources  of  irradiation. 

The  model  described  here  is  a  first  approximation  to  a  final  model 
that  will  be  used  to  predict  the  characteristics  of  naturally  occur¬ 
ring  radiation  incident  on  a  target.  The  work  described  here  represents 
the  status  of  this  portion  of  the  Ground  Target  Signature  program  six 
months  after  the  inception  of  the  program. 

II.  PROBLEM  DESCRIPTION 

The  effort  described  in  this  report  is  the  first  step  in  predict¬ 
ing  the  characteristics  of  radiation  incident  on  a  ground  target.  In 
order  to  simplify  the  process  at  this  stage,  several  limiting  features 
were  incorporated  into  the  model.  First,  the  wavelengths  being  con¬ 
sidered  in  this  report  extend  from  0.4  to  1.7  microns  in  0.0S  micron 
steps.  This  region  was  selected  since  it  is  of  greatest  concern  at 
the  present  time,  and  since  limited  validation  could  be  performed  in 
the  visual  spectral  region  with  equipment  already  on  hand.  Second, 
the  situations  being  considered  in  this  report  limit  the  input  to  those 
from  the  sun  and  the  cloudless  daytime  sky.  The  elimination  of  clouds 
from  the  problem  precludes  the  immediate  need  to  mathematically  des¬ 
cribe  a  cloud  and  simplifies  correction  and  validation  of  the  pre¬ 
liminary  model.  Third,  radiation  incident  on  the  target  after  re¬ 
flection  from  the  ground  is  not  being  considered  at  this  time  since 
its  effect  is  much  less  than  that  of  the  sun  or  sky.  Also,  except 
for  very  simple  terrestrial  configurations,  the  amount  and  distri¬ 
bution  of  radiation  from  this  source  would  be  much  more  difficult  to 
model.  Once  a  basic  model  has  been  generated  and  validated,  these 
perturbing  influeces  can  be  added  separately  to  create  a  more  general 
model  that  can  be  applied  for  any  set  of  precisely  defined  atmos¬ 
pheric  and  terrestrial  conditions. 

■K 

References  are  listed  on  page  32. 
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The  model  described  in  this  report  predicts  the  spectral  irra- 
2 

diance  (watts/un  in  a  0.C5  micron  band)  reaching  the  ground  from  the 
sun  and  the  clear  daytime  sky.  The  solar  input  is  the  solar  irradiance 
outside  the  atmosphere  multipled  by  the  atmospheric  transmission  to 
the  ground.  The  spectral  irradiance  produced  by  a  defined  segment 
of  sky  is  determined  through  the  format  indicated  in  Figure  1. 

A  direct  solar  ray  of  irradiance  enters  the  atmosphere,  and  a 
fraction  T^  of  this  irradiance  reaches  a  volume  element  within  the 
field  of  view  of  a  ground-based  sensor.  A  certain  fraction  of  the 
radiation  will  be  scattered  within  this  element.  The  amount  of  energy 
scattered  within  the  element  and  the  manner  in  which  the  scattered 
radiation  is  distributed  in  space  are  determined  from  the  total  and 
directional  volume  scattering  functions.  The  radiation  scattered 
within  the  cone  il  is  then  attenuated  along  the  path  to  the  ground. 

The  transmission  along  this  leg  is 

The  total  power  ai living  at  the  sensor,  assumed  for  convenience 

to  have  an  area  of  one  square  centimeter,  in  the  wavelength  band  AX  is 

« 

the  sum  of  the  contribution  from  all  volume  elements  within  the  field 
of  view.  The  result  is 

H  =  Hs  J  B(6) OTjT2dV  ,  (1) 

r=o 

where  r  is  the  distance  from  the  sensor  to  the  scattering  volume.  Hg 
can  be  removed  from  the  integral  since  it  has  no  r  dependence.  This 
equation  is  strictly  valid  only  for  single  scattering.  Multiple 
scattering  occurs  when  the  sky,  in  addition  to  the  sun,  is  a  source  of 
radiation  which  can  be  scattered  into  the  field  of  view.  It  is  of 
particular  importance  at  short  wavelengths  where  the  scattering  ef¬ 
ficiency  is  high.  The  effect  of  multiple  scattering  will  be  addressed 
later  in  the  program.  The  next  section  will  be  devoted  to  determining 
analytical  expressions  in  terms  of  r  for  the  various  terms  in  equation 
(1)  based  on  physical;  geometric,  and  environmental  parameters. 
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III.  EVALUATION  OF  MODEL  INPUT  VARIABLES 

A.  Evaluation  of  H 
_ s 

is  determined  by  an  integration  of  the  Planck  radiation  equation 
and  can  be  expressed  as 

(2) 

2  .  2 

Hs (watts/cm  in  AX)  =  ir(sr)  x  N(watts/cm  -sr)  x  D(fraction  of  N  in  AX). 


The  quantity  ttN  is  a  direct  function  of  solar  temperature  (5750  degrees 

K) ,  solar  radius  (6.956  x  10*°  cm),  and  solar  distance  (1.495  x  lO1^  cm) 

2  3  2 

The  emitted  solar  irradiance  is  taken  as  6.211  x  10  watts/cm  .  The 

2 

use  of  this  value  results  in  uN=0.1344  watts/cm  .  This  quantity  is  the 

2 

solar  constant,  and  is  usually  given  in  units  of  calories/cm  -min. 

t  4 

Other  references  give  its  value  as  0.1350,  0.1393  i  0.0028,  and 

5  2 

0.1392  watts/cm  .  For  calculations  in  this  report  we  will  use  the 

2 

value  7rN  =  0.139  watts/cm  . 


The  fraction  of  this  irradiance  in  various  wavelength  oands  is 
given  by 


X+AX  X  /•» 

° '[  A“  -  Adx]  /S 


(3) 


-2  -1  .  -1 

where  (watts  cm  sr  micron  )  is  the  spectral  radiance.  The 
integral  ratios  given  in  equation  (3)  are  tabulated  in  Reference  2. 

By  assigning  values  to  X  (0.40,  0.45,  0.50,  etc.,  microns)  and  AX 
(0.05  microns),  we  can  determine  H  ,  as  tabulated  in  Table  I.  These 
results  are  also  shown  in  Figure  2  to  provide  a  more  graphic  description 
of  the  wavelength  dependence  of  the  function. 


B.  Evaluation  of  fldV 

Both  Q  and  dV  can  be  determined  from  geometric  parameters.  The 

value  of  0  is  approximated  by  A/r^,  where  A  is  the  area  of  the  sensor. 

2  -2 

Since  we  are  assuming  a  sensor  area  of  1  era  ,  Q  =  r  .  The  scattering 

2  2  d 

volume  is  given  by  dV  =  nr  tan  *  y  dr,  where  e  is  the  sensor  field  of 
view.  The  product  OdV  is  then  given  by  it  tan^  y  dr.  The  final  express! 
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is 

fldV  =  10“5  *  tan2  |  dr,  (4) 

where  dr  is  now  expressed  in  kilometers. 

C.  Evaluation  of  B(6) 

The  directional  scattering  efficiency  B(6)  is  a  measure  of  the 

efficiency  with  which  the  incident  irradiance  is  scattered  into  the 

3 

sensor  direction  by  1  cm  of  dV.  It  is  a  function  of  wavelength, 
the  scattering  angle,  and  r,  since  the  number  of  scattering  particles 
per  cubic  centimeter  decreases  with  altitude.  Two  expressions  for  B(6) 
are  required,  since  radiation  can  be  scattered  either  by  Rayleigh  or 
Mie  particles. 

The  first  step  in  determining  B(6)  is  calculating  6  as  a  function 
of  the  geometric  parameters.  The  expression  which  results  is 

cos  6  *  cos  6  cos  0.  cos  <f>  +  sin  6.  sin  0_,  (5) 

S  L  L  S 

where  0g  and  are  the  elevation  angles  of  the  sun  and  the  sensor 
field  of  view  respectively,  and  <p  is  the  azimuthal  angle  between  them. 
These  angles  are  specified  as  input  parameters.  The  elevation  3nd 
azimuthal  angles  of  the  sun  can  be  found  in  standard  references, 
equation  (5)  is  plotted  in  Figures  3  to  6  for  various  values  of  0c, 

0L,  and  4>. 

Raleigh  scattering,  which  is  a  special  case  of  Mie  scattering, 

occurs  when  electromagnetic  radiation  impinges  on  particles  that  are 

significantly  smaller  than  the  wavelength  of  the  radiation.  In  this 

problem  the  Rayleigh  particles  of  importance  are  air  molecules.  For 

unpolarized  incident  radiation,  the  equation  describing  the  Rayleigh 
.  7 

scattering  is 


(6) 


2  2 

2 n  (n-1)  ,,  2„^ 

B(6)p  ~  - 4 — —  (1  ►  cog  6) 

R  NX4 


where  n  is  the  index  of  refraction  of  air  and  N  is  the  number  of 

3 

molecules  per  or.  .  The  wavelength  dependent  portion  of  this  expression 
fr-112 

is  ,  which  is  given  for  standard  temperature  and  pressure  a?  a 

A 

function  of  wavelength  in  Table  1.  This  information  was  extracted 

8  19 

from  Reference  7  which  cited  Kuiper.  N  can  be  shown  to  be  2.67  x  10 

3 

molecules  per  cm  at  zero  degrees  centigrade  and  pressure  of  1013 
millibars. 


The  altitude  dependence  of  8(6)  is  determined  from  ideal  ga.s  laws 
assuming  an  isothermal  atmosphere.  Both  n-1  and  N  decrease  with 

V/g  y 

altitude  as  e~n/  ,  where  h  is  the  altitude  in  kilometers.  The  scale 
height  of  the  Rayleigh  atmosphere  is  8  km,  which  is  the  altitude  at 
which  the  density  of  the  atmosphere  drops  to  one  third  its  sea  level 
value.  It  is  also  the  total  height  the  Rayleigh  atmosphere  would  have 
if  it  were  all  at  sea  level  density.  Our  final  expression  is  therefore 


2iTxl05  (n  -l)2  „  ,  .  WQ  . 

B(6)  =  - -  (1  +  cos  5)  e  L  km  sr  .  (7) 

K  N  X4 
o 


The  subscript  R  indicates  that  it  is  the  Rayleigh  component  of  B(6); 

the  subscript  c  indicates  sea  level  values  of  these  quantities.  The 

2 

only  6  dependence  arises  through  the  term  1  +  cos  5.  Therefore, there 
is  a  two-to-one  difference  between  the  maximum  and  minimum  directional 
scattering  efficiencies  in  a  Rayleigh  medium.  This  will  be  seen  to 
contrast  sharply  with  the  properties  of  Mie  scattering  that  will  be 
discussed  next. 
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Mie  scattering  occurs  whea  ladiation  impinges  on  particles  whose 
radius  is  greater  than  or  equivalent  to  the  wavelenrth  of  the  radiation. 
The  total  Mie  scattering  coefficient  is  given  oy' 


oc 


(8) 


In  this  equation  it  is  assumed  that  the  droplets  are  spherical.  This 

is  a  good  assumption  for  droplets  comprising  haze  or  fog.  nR  is  the 

number  of  iroplets  per  cubic  centimeter  with  radii  between  R  ana  R+dR, 

and  K„  is  the  scattering  area  ratio  which  is  a  single-valuei  function 
8  9 

of  2ttR/a.  Kr,  which  was  first  given  by  Houghton  and  Chaliccr,  is 

shown  without  fine  structure  in  Figure  7.  Physically,  is  the  ratio 

between  the  optical  scattering  cross  section  and  the  geometric  cross 

section  of  a  particle.  Dessens1®  has  estimated  that  the  mean  calue 

of  R  in  continental  hazes  is  0.4  micron  for  a  high  (*v  78  percent) 

relative  humidity.  If  we  assume  that  all  droplets  have  this  radius 

and  that  the  total  Mie  coefficient  can  be  converted  to  a  directional 

Mie  coefficient  by  introducing  a  function  of  6,  equation  (8)  becomes: 

B(5)m  ^  S.,o:  *  10~4  n  f (6)  km"1  sr"1  .  (S) 


Kg  is  given  in  Table  I  as  a  function  of  X  for  R  =  0.4  micron. 


tj  can  be  inferred  from  the  horizontal  visibility  V  through  the 
following  procedure.  Hie  total  scattering  coefficient  for  a  pre¬ 
dominantly  scattering  atmosphere  is  given  by  8  =  3.912/V  11 ,  where  8 

is  the  ;um  of  the  Rayleigh  and  Mie  total  scattering  coefficients  in 

-4 

the  visual  region.  The  total  Mie  coefficient  BT  M  Q  =  5.03  y  20 
n  KR  km 

equation  (6)  over  the  complete  sphere  resulting  in 

2  r  .  >  2  .o  3. 


The  total  Rayleigh  coefficient  can  be  found  by  integrating 


,R,0  =  f 


2^)1  (1  +  cos2  6)  d,  . 


km 


-1 


(10) 


4ir 

steradians 


NX 


3NX 
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This  coefficient  is  listed  in  Table  I  as  a  function  of  X  for  the 
visual  portion  of  the  spectrum. 

The  value  °f  ^  q  can  now  expressed  as 

6t,m,o=°E  gjs-^r--  \R,o)  (11) 

X=0 . 4  ’ 

for  the  visual  portion  of  the  spectrum.  The  function  is  a  weighting 
factor  representing  the  wavelength  response  of  the  human  eye.  The  eye 
response  is  shown  in  Figure  8  plotted  from  data  given  in  Reference  2. 
The  curve  was  normalized  to  unit  area,  and  the  fractional  areas  in  0.05 
micron  steps  were  determined.  These  values  of  g  are  given  in  Table  I. 

A 

A  visibility  can  be  assigned  resulting  in  a  value  of  8^.  ^  ^  from 
equation  (11).  When  combined  with  KR  in  equation  (8),  n  =  n(V)  can 
be  determined  with  the  results  plotted  in  Figure  9. 

This  procedure  can  now  be  applied  for  other  wavelengths.  A  value 

of  V  is  assigned  to  determine  n  and  a  wavelength  is  assigned  to  deter- 

-4-1 

mine  K„.  Then  $„  ,,  is  5.03  x  10  n  L  Kd  .  The  results  are  shown 
R  T,M  R 

in  Figure  10. 

The  angular  distribution  of  the  Mie  scattered  radiation  can  be 

determined  as  follows.  The  total  power  in  AX  scattered  from  one  cubic 

centimeter  of  Mie  particles  is  Kg  8^  ^  x  10  5  watts.  The  angular 

distribution  of  this  power  is  given  by  curves  found  in  several  places 

1 1  12 

in  the  technical  literature,-''*  The  form  used  in  this  report  u as 
taken  from  Dcirmendjian. *  These  data  were  plotted  as  f(6)  in 
Figure  il  after  normalization  such  that 

j  f(5)dn  a  1  . 

4tt 

steradians 
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The  variation  of  8p  ^  *  8p  ^  with  a  altitude  is  indicated  by 
11  14 

Middleton  ,  citing  Waldram.  The  curves  are  shown  in  Figure  12.  The 
letters  assigned  to  some  of  these  curves  will  be  discussed  in  the  next 
section.  It  can  be  seen  that  there  are  large  differences  between  the 
various  curves.  On  the  average,  the  Mie  volume  scattering  coefficient 

varies  with  altitude  as  e"^/4.1  for  h  >  one  kilometer;  however,  it 

should  be  noted  that  for  high  ground  level  values  of  Bp  M  the  rate  of 

decrease  of  £„  w  with  altitude  is  much  more  rapid  at  the  lower  altitude 
1 

than  at  the  higher  altitudes.  The  model  used  to  approximate  these 
curves  will  be  discussed  in  more  detail  in  the  next  section,  but 
briefly  the  situation  is  as  follows.  If  Bp  is  small  [V  >  G(A}],  the 

altitude  variation  goes  h/4.1km.  j  f  v  <  G(A),  the 

curves  are  approximated  by  two  line  segments  such  that  the  scale  height 
is  4.1  kilometers  for  h  >  1  km  and  a  kilometers  for  h  <  1  km.  We  will 
now  determine  values  for  G(A)  and  a. 


G(A)  can  be  determined  by  specifying  the  value  of  Bp  M  at  some 
altitude.  This  specification  was  made  somewhat  arbitrarily  by  an 

assumption  that  Bp  M  yisuai  region  <  0*1  at  331  altitude  of  one  km. 

The  curve  passing  through  the  Bp  ^  =  0.1  km-1,  h  =  1  km  point,  with  a 

scale  height  of  4.1  km,reaches  the  zero  altitude  line  at  ^  =  0.128  km  1 

The  visibility  represented  by  regi(m  =  0.128  k.'1  is 

designated  G(A)  and  can  be  found  by  replotting  Figure  10.  The  result 
is  shown  in  Figure  13. 


The  scale  height  of  the  lower  atmosphere,  a,  can  be  determined 
by  finding  the  equation  of  the  exponential  curve  connecting  the  points 

0.7 

ftp  M  =  E  5.03  x  10"4  nKj^,  h  =  0)  and  (8T  M  =  0.1  km'1,  h  =  1  km). 
A=0.4 


The  result  is  a  =  -1/ln  (0.1/Bp  H  ,  where  @T  M  0 


is  the  Mie  volume 
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scattering  coefficient  at  ground  level  shown  in  Figure  10.  After 
substituting  h  =  r  sin  0^,  the  following  expressions  <  an  be  applied 
directly. 


If  V  >  GM,  BTjH  »  6TjM>0  e-«r  sin  V/4-1  0  <  r 


(12) 


If  V  <  G(A), 


ST,M  *  6T,M,0  e  T  Si"  °L  l"(°'1/liT,l  ,0>  0  <  r  <  1/sin  eL, 


^0T  M  =  0.128  e‘fr  sin  Q0/4A  l/si  .  0L  <  r  <  -  , 


where  &]-  ^  q  =  5.03  x  10  4  n  K^. 
coefficient  is  then  given  by 

B(6)m  =  5.03  x  10 


The  directional  Mie  scattering 
'4  n  ^  f(5)  B(r)  , 


(13) 


where 


if  V  >  G(A), 


if  V<  G(A), 


B(r)  =  e-(rsi"  8l)/4-i 

B(r)  =  e  r  sin  ®L  ta(0-1/6T,H,0) 

B(r)  =  e‘tr  sin  eL)/4-1 


0  <  r  <  00 
0  <  r  <  1/sin  0L 
1/sin  <  r  <  00 


(14) 


D.  Evaluation  of 

The  transmission  over  the  two-legged  patli  (Figure  1)  is  the  product 
of  foi’"  individual  transmissions.  The  transmissions  through  the 
Rayleigh  and  Mie  scattering  components  have  already  been  mentioned. 
There  is,  however,  also  absorption  by  water  vapor  and  carbon  dioxide, 
which  will  be  discussed  using  the  geometry  shown  in  Figure  14. 

The  horizontal  transmission  by  the  Rayleigh  component  of  the  atmo¬ 
sphere  is  Tr»  =  exp  (-r*  8^.  R).  However,  for  slant  paths,  8T  R  is  a 
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function  of  r'  so  we  must  evaluate  the  egression  Tr  =  exp  I-  |  B(r')dr' 

*-  o 

To  determine  the  form  of  we  must  integrate  equation  (7 ,  with 

respect  to  dfi  over  all  directions.  The  result  is 

,,  „  32,l(a-li:x.ipL  e-Cr'  sin  ep/8  k„,-l  (ls, 

1,K  3NX4 

When  this  value  is  substituted  into  the  transmission  equation  using 
the  geometry  of  Figure  14,  the  result  is 


rr.  eT,,,,o  e'(r'  sin  ;)/S  dr'J 

Ur  Km*- <r ^ 4 * sin es)/8 j i - 

1  U  A 


,  (16) 


When  the  integrations  are  performed  and  the  r  dependence  is  factored  out, 
the  resulting  expression  is 


where  R  0  is  given  by  equation  (10) . 

The  absorbing  components  of  the  atmosphere  behave  in  a  different 
manner  from  the  scattering  components.  For  convenience,  we  will  define 
a  distance  which  contains  the  same  amoiait  of  absorbing  material  as  the 
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total  light  path.  This  distance  will  be  called  the  optical  slant  range. 

The  carbon  dioxide  optical  slant  range  will  be  given  the  symbol  p. 

It  will  be  assumed  here  that  since  carbon  dioxide  is  a  permanent  atmos¬ 
pheric  gas  its  concentration  will  decrease  with  altitude  at  the  same 
rate  as  the  Rayleigh  component.  Thus,  using  the  geometry  of  Figure  14, 

T  00 

p  =  (  e-(r'  Sin  0L)/S  dr1  +  |  e'(r  Sin  6L  +  R  Sin  6S)/8  dR  .  (18) 

o  o 

When  the  integrations  are  performed  and  the  r  dependence  is  factored  out, 
the  resulting  expression  is 


P  = 


sin  8. 


-  SC-U-j - )  e'(rsin  eL)/Skm  . 


sin  0g  sin  eL 


This  equation  is  plotted  in  Figures  15,  16,  and  17.  Several  features 
of  these  graphs  are  readily  apparent,  as  listed  below. 

(1)  When  6L  =  8g,  p  is  dependent  on  r. 

(2)  When  8^  <  8g,  p  increases  with  increasing  r. 

(3)  When  QL  >  g,  p  decreases  with  increasing  r. 

8  8 

(41  p  =  — — —  or  — — —  ,  whichever  is  greater. 
v  •  max  sin  8g  sin  SL  *  6 

(5)  pjnin  =  8  kilometers,  the  CX>2  atmospheric  optical  depth. 

It  has  been  demonstrated*8'*6  that  carbon  dioxide  absorption  can 
be  described  by  an  error  function  absorption  law,  namely 

^  2 

T  =  1  -  erf  K(X)  /p,  where  erf  Z  *  ~  [  e  t  dt  .  (2C 

c  /n  i 

o 

K(A)  as  used  in  this  study  was  determined  by  plotting  transmissicn 
data  given  in  Hack  forth,*5  integrating  under  the  curve  in  0.£b  micron 
bands,  and  c*>n**arting  the  results  to  obtain  a  value  of  K(A)  for  that 
band. 


20 


— .  ... 


The  error  function  can  be  expanded  in  an  infinite  series  as 


erf  Z  = 


4tstirin-+  ^ 


The  maximum  value  of  K(X)  in  the  0.4  to  1.7  micron  interval  is 
0.007S.  The  maximum  value  of  p  will  be  taken  as  8/sin  10°  =  46.07. 

Thus  Zfflax  =  0,0075  A6.07  =  G.050P.  Therefore, 


erf  Z_ 


2x0.0509 


0,00259 

3 


0.0000067 

10 


...]  .  (22) 


Obviously.,  the  second  and  succeeding  terms  can  be  ignored;  thus  our  final 
expression  for  carbon  dioxide  transmission  is 


Tc  =  1  -  =-  K(X)  1  -  a(X)/p  ,  (23) 

Jr 

where  p  is  given  by  equation  (19).  The  values  of  a(A)  are  given  in 
Table  I. 

Absorption  by  water  vapor  follows  a  similar  law,  namely,15 

TH  =  1  -  erf  M  ,  (24) 

where  8(A)  is  the  water  vapor  error  function  absorption  coefficient 
and  X  is  the  number  of  centimeters  of  precipitable  water  vapor  in  the 
light  path.  The  distribution  of  water  vapor  with  altitude  can  be 
found  on  a  graph3  relating  the  altitude  and  the  mixing  ratio  (grams  of 
water  per  kilogram,  of  dry  air)  for  a  ground  level  humidity  of  7.35 
grains  of  water  per  cubic  meter  of  dry  air.  This  curve  was  altered  to 
represent  W*  (cm  of  precipitable  water  vapor  per  km  of  horizontal 
path)  as  a  function  of  altitude.  The  results  are  shown  in  Figure  18. 
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It  can  be  seen  that  on  the  logarithmic  plot  the  low  altitude 
portion  of  the  curve  can  be  approximated  by  a  straight  line  equivalent 
to 


W*  =  W*  e 
o 


-h/2 


(25) 


where  W*,  the  ground  level  value  of  W*,  is  a  function  of  temperature 
and  relative  humidity  as  shown  in  Figure  19. 

W  can  be  determined  using  the  geometry  of  Figure  14.  Thus, 

r 

W  =  J  W*  e'ir'  sin  V'/2  dr’ 

o  (26) 

00 

♦  |  W*  e’(r  Sin  eL  +  R  Sin  eS)/2  dR  . 
o 

When  the  integrations  are  performed  and  the  r  dependence  is  removed,  the 
result  is 


2W 


W  = 


sin  6, 


2K  'ii H 


1 


sin  0, 


e-(r  sin  6l)/2 


(27) 


This  equation  is  plotted  in  Figures  20,  21  and  22.  Again,  several 
features  of  these  graphs  are  readily  apparent: 


(1) 

When 

«L='S 

(2) 

When 

8L  ‘  9S 

(3) 

When 

8.  >  8C 
t-t  o 

2W* 

(4) 

W 

o 

max 

sin  8, 

i 

(5) 

W  . 
mm 

=  2W£  , 

W  is  independent  of  range. 

W  increases  with  increasing  range. 

W  decreases  with  increasing  range. 

2W* 

‘  or  •  q"  ,  whichever  is  greater, 
sin 

the  water  vapor  atmospheric  optical  depth. 


The  absorption  coefficient  was  determined  by  integrating  trans¬ 
mission  data  given  in  Hackforth.^  The  maximum  value  was  found  to  be 

1.264.  The  maximum  value  of  W  will  be  taken  as  -----  — -  =  116.3. 

sin  10° 

Thus  Z  =  0.632  Al6. 3tt  =  12.1.  Since  the  error  function  series 
max 

expansion  converges  very  slowly  for  this  value  of  Z,  a  rational  approxi- 
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mation  to  erf  Z  was  found  *  for  ease  in  programming  the  equation. 

This  approximation  states: 

T  _  1  -  7  f  0.34802  0.09588 

H  -  1  '  erf  Z  -  L(i+o.47047Z)  “  (1+0. 47047Z)2 

(28) 

0.74786  1  -Z2  . 

+  - -  e  +  e(z)  > 

(1+0.47047Z)JJ 


where  |e(Z)|  <  2.5x10  ^  and  Z  =  -  ^  /irW  =  y(A)  JH.  y(A)  is  given  in 
Table  I. 


The  transmission  through  the  Mie  scattering  component  agr.in 
requires  separation  on  the  basis  of  the  ground  level  visibility.  If 
V  >  G(A),  we  car  use  equation  (32)  and  the  geometry  of  Figure  14.  The 
result  is 


tm  “  exP  - 


Si"  V/4a  *'! 


[29) 


exp  -II  6i,m,o  e‘tr  sln  'L  *  R  Sin  8s,/4'1  * 


This  expression  reduces  to 
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If  V  <  G(X),  we  must  distinguish  between  the  cases  where  r  <  1/sin  0^ 
and  r  >  1/sin  0^.  The  geometries  associated  with  these  cases  are  shown 
in  Figures  23  and  24.  The  values  of  6^.  ^  for  each  segment  are  also 
indicated  on  these  Figures.  When  r  <  1/sin  0^  (Figure  23), 


tm  ' 


j-  I  0.12 
1  1 


28e-(Rsin  V74-1  dR  >  x 


sin  6, 


exp 


.  sin  0C 

■I 

'  r  sin  0. 


t  eR  sin  ®S  lnt°'1/6>  dsj  x 


(31) 


sin  0, 


exp 


I- !  ■  -  “ 


sin  0^  In (0.1/&)  ^ 


i  • 


where  0T  M  Q  is  given  the  symbol  0.  This  expression  reduces  to 


Sj 


.  exp 

^  j  sin  0L 


0.1/ln (0.1/3) 

sin  9g 


♦  0.41 


3 

f  1  1  1 

In (0.1/3) 

sin  6g  sin  0^ 

sin  0^ 


0,  In (0.1/ 


When  r  > 1/sin  0L  (Figure  24) , 


tm  ■  exP 


!• ! 


0. 128  e_{-R  Sin  V/4 


r  sin  0j 
sin  0,, 


A  dR  j 


sin  6, 


exp 


(9X12 
■  ! 


(  (  n  l'yo  ” (r  sin  0,)/4.1  ,  •) 

exp  <  -  0.128  e  v  L  dr  > 


(33) 


sin  §. 


r'  sin  0.  InfO 
e  L 


(o.im  dr,j 
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This  expression  reduces  to 


t  = 


M 


=  exp 


exp 


(8-0.  l)/ln(0. l/g~)  -  0.41 


sm  6, 


i 

1 


(34) 


| -0.51  |^sin  -  sin  e^J 


•(  sin  6L)/4.1 


These  expressions  can  be  checked  in  part  by  evaluation  at  the 
limits  of  applicability.  If  equations  (32)  and  (34)  are  evaluated  at 
the  point  r  =  1/sin  0^,  they  are  identical  expressions,  namely. 


M 


exp 


(8-0.1)/ In (0.1/8) 
sin  eL 


0.41 

sin  es  ;  ‘ 


(35) 


If  r  =  1/sin  0L  is  substituted  into  equation  (30)  and  equations  (30) 
and  (35)  are  solved  at  the  limit  V  »  G(A)  (8  =  0.128),  identical 
expressions  result,  namely. 


TM=eXP 


-0.11 

sin  8, 


0.41 
sin  0- 


(36) 


One  other  check  is  available.  If  we  let  0^  =  0g  =  90  degrees  in 
equation  (36),  we  should  arrive  at  the  same  solution  as  the  horizontal 
transmission  solution  of  TM  =  exp  (-R8),  where  R  =  4.1  km  (the  Mie 
scattering  scale  height  for  V  =  G(A)  and  8  =  0.128).  This  also  is 
found  to  check,  resulting  in  TM  =  0.59. 
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IV.  UTILIZATION  OF  THE  MODEL 

As  shown  in  Section  II,  the  problem  can  be  reduced  to  the  eval¬ 
uation  of  an  integral  containing  several  r-dependent  variables  over 
all  values  of  r.  The  integral  in  its  most  usable  form  is  obtained  by 
rewriting  equation  (1)  as 

r=* 

H  =  Hg  |  B(6)0TrTmTcTh  eft'  .  (37) 

r=o 


Each  of  these  variables  has  beer,  evaluated  as  a  function  of  r  and  X 
earlier  in  this  report.  For  convenience,  the  location  of  each  will 
now  be  summarized.  Hg  is  given  in  Table  I.  The  geometric  variable 
fldV  is  given  by  equation  (4).  B(6)  is  the  sum  of  the  components  from 

Rayleigh  and  Mie  scattering.  The  Rayleigh  component  is  described  by 
equation  (7)  where  Nq  =  2.67  x  1019,  (nQ-l)2/A4  is  given  in  Table  I, 
and  6  can  be  determined  from  equation  (5) .  The  Mie  component  is  de¬ 
scribed  by  equation  (13)  where  n  is  given  in  Figure  9,  KR  is  given  in 
Table  I,  f(6)  is  shown  in  Figure  11,  and  B(r)  is  given  by  equations (14) . 
G(A),  the  visibility  parameter,  is  given  in  Figure  13. 

Hie  transmission  requires  four  terms,  since  scattering  by  Rayleigh 
and  Mie  particles  and  absorption  by  carbon  dioxide  and  water  vapor 
must  be  taken  into  account.  TR  is  described  by  equation  (171  where 
r  o  is  given  in  equation  (10).  TM  is  described  by  equations  (30), 
(32)  and  (34)  with  G(A)  again  given  by  Figure  13.  3-r  u  n  defined 

as  5.03  x  10~^  riKR(X).  Tc  is  described  by  equation  (20)  with  a(A) 
given  in  Table  I  and  p  given  by  equation  (19).  is  given  by  equation 
(28)  where  Z  =  y(A)  y(A)  is  given  in  Table  I,  W  is  given  in 
equation  (27),  and  W*  can  be  determined  from  Figure  19. 

At  upper  limit  of  infinity  cannot  be  handled  by  the  computer; 
therefore,  an  upper  limit  was  assigned  based  on  the  following  procedure. 
The  most  slowly  converging  term  in  the  integral  is  the  Rayleigh 
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scattering  term  having  a  scale  height  of  eight  kilometers.  We  can 
therefore  determine  the  value  of  b  for  which 

b 


sin  0 


8 


I6'" 


L  [  _-r  sin  0l/8  dr  =  Q>999 


and  be  assured  that  substitution  of  this  value  for  the  upper  limit  in 
equation  (37)  will  result  in  an  error  in  the  value  of  the  integral  of 
less  than  0.1  percent.  This  value  of  b  is  found  to  be  56/sin  0^. 

These  terms  can  be  combined  and  simplified,  resulting  in  an 
equation  of  the  form 

Q 

H  =  A  j  {DeFr  +  GeHr}  {exp  [i  *  JeHr]  f  jexp  [l(eFr]j  x 


L  I  JuJ*  fr  0.34802  0.09588 

r  '  L  "fe  -  J'F/  | [^(1+0.470472)  “  (1+0.47047Z)2  (38:) 


0.74786  1  -Z2\  . 

U+G.47047Z)Jj X  e  f  d  ’ 


where  A  =  tan 


‘-9. 

f  exp  - 


-9.9  x  10'12  (n  -1)2/A4] 

~&nr - }  D  - 2-3  *  10  x 


(l+cos26)  (nQ-l)2/A4,  F  *  sin  0j/8,  G  *  1.58  x  10'8  n^ffS), 
K  =  9.9  x  10~12  (CSC  0S  -  CSC  eL)  (n  -1)2/X4,  L  *  a(A), 


M  =  8 (esc  Sj.  -  esc  0^),  and 


Z  =  y(A) 


/  2W* 

VsITi -+2w:  Ccsc 


a  „  .  -r  sin  8,  /2 

0g  -  esc  0^)  e  L 


27 


£i*2CSX&»i* 


The  remaining  inputs  are  functions  of  G(a),  If  V  >  G(>J,  only 
one  integral  is  required  and  B  =  0,  C  =  S6/sin  6^,  H  =  -  sin  6^/4. 1, 

I  =  --4.13/sin  0L,  and  J  =  -4.1g(csc  Og  -  esc  0^) .  If  V  <  G(A),  the 
required  solution  is  the  sum  of  two  integrals.  The  first  is  rep¬ 
resented  by  B  =  0,  C  =  1/sin  0L,  H  =  sin  eL  ln(Q.l/g), 

I  =  8/lnCO.l/g)  _  0.1/ln(0.i/g)  +  0.41  on 
sin  6^  sin  0g  * 

J  =  g(csc  eg  -  esc  0^)/ln(O. 1/g) .  The  second  integral  is  represented 


by  B  =  1/sin  6. ,  C  =  56/sin  0,  ,  H  =  -sin  0./4.1, 


T  (g-0. 1) / In (0. 1/6) -0.41  ,  T  „ 

I  =  n - L - »  and  J=  -0.51  (esc  0C  -  esc  0t). 

sin  o  s  1 


appears  in  these  expressions  is  6^  M  q  given  by  5.03  x  10 


The  g  which 
nKR. 


In  order  to  provide  a  limited  validation  of  part  of  this  model,  a 
comparison  was  made  between  the  transmission  predicted  by  the  nxyiei. 
from  the  ground  level  visibility  and  that  determined  £rc;a  actual  values 
of  the  altitude  dependence  of  the  visual  spectrum  volume  scattering 
coefficient.  The  Mie  transmission  predicted  by  the  model  xs  given  by 


M 


The  value  of  B..(r)  for  this  expression  is 

"I 


o 

shown  in  Figure  25  as  a  function  of  the  ground  level  visibility.  ^her 
V  <  G(A),  two  curves  are  required.  As  V  approaches  G(A),  the  slope  of 
the  lower  segment  approaches  that  of  the  upper  segment.  At  V  =  G(A) 
only  one  curve  is  required,  and  for  further  increase  in  V  the  entire 
curve  is  translated  while  keeping  the  same  slope.  The  Mie  transmission 
along  one  of  these  curves  to  an  altitude  of  nine  kilometers  is  found 
by  letting  0.  equal  90  degrees  and  integrating  equation  (13)  from 
r  =  0  to  r  =  9  km.  6^  ^  ^  is  given  by  equstion  (11).  The  Rayleigh 
transmission  can  be  found  by  letting  9^  equal  90  degrees  and  integrat¬ 
ing  equation  (15)  from  r  o  0  to  ;  =  9  km.  The  product  of  the  Rayleigh 
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ana  Hie  transmissions  represents  the  total  transmission  ever  the  verti¬ 
cal  path  to  an  altitude  of  nine  kilometers  as  predicted  by  the  model. 

This  procedure  was  applied  to  the  curves  labeled  (a)  through  {£} 
in  Figure  12,  The  individual  curves  were  also  approximated  by  a  large 
020)  number  of  straight  line  segments  r  the  transmission  over  each 
segment  was  determined.,  said  the  product  of  the  transmissions  was  found 
for  each  curve.  The  transmissions  determined  by  the  two  methods  were 
compared  and  errors  of  3.9,  2.4,  1.2,  27.5,  0.8,  and  10.5  percent  were 
found  for  curves  (a)  through  (f)  respectively.  The  large  error  on 
curve  (d)  resulted  from  the  increase  in  &T  ^  with  altitude,  which  is 
not  a  common  occurrence  on  a  clear  day. 

V.  CONCLUSIONS 

An  atmospheric  radiative  transfer  model  based  on  single  scattering 
has  been  described  in  this  report.  At  this  stage  it  is  not  known  h.  * 
useful  it  is;  however,  certain  features  are  clear.  First,  the  geometry 
shown  in  Figure  1  is  generally  applicable  to  problems  of  this  kind, 
provided  multiple  scattering  can  be  ignored.  Second,  equation  (1) 
fully  describes  the  geometry  in  the  figure,  although  defining  the 
terms  in  the  integral  for  non- ideal  conditions  will  not  generally  be 
so  straightforward  as  described  hers.  Third,  it  is  desirable  to  use  a 
simple  model  of  this  nature,  sines  at  the  time  a  data  point  is  desired 
it  is  unlikely  that  detailed  information  about  the  scattering  and  ab¬ 
sorbing  properties  of  the  upper  atmosphere  will  be  known.  Thus,  some 
measure  of  predictability  based  on  measurements  that  can  be  made  at 
ground  level  is  required. 

In  the  course  of  this  program,  several  sources  of  error  have  been 
found.  The  errors  tend  to  occur  in  three  general  classes.  The  first 
class  of  errors  results  from  simplifications  and  approximations  used 
to  keep  the  model  from  being  overly  complicated.  One  such  error 
results  from  the  atmospheric  transmission  model  being  based  on  a  flat 
earth  assumption.  The  validity  of  this  assumption  is  a  strong  function 


29 


of  9,  and  9  .  An  analysis  was  made  of  this  error  source,  and  it  was 
found  that  the  error  is  less  than  one  percent,  provided  that  both  eL 
and  are  greater  than  five  degrees.  Another  error  in  this  class 
results  from  the  manner  in  which  the  Mie  scattering  coefficient  is 
derived  from  the  visibility.  Tne  Mie  scattering  coefficient  used  in 
the  model  is  based  on  the  number  of  droplets  per  cxbic  centimeter  in 
the  volume  element.  It  is  assumed,  for  convenience,  that  all  droplets 
have  the  same  radii,  corresponding  to  a  mean  radius  found  in  continental 
hazes.  In  reality,  however,  there  is  a  spread  in  radius.  The  effect 
of  this  error  is  currently  being  analyzed. 

The  second  class  of  error  results  from  inadequate  knowledge  of 
some  of  the  basic  input  parameters.  Some  of  these  uncertainties  in¬ 
clude  the  manner  in  which  the  concentrations  of  water  vapor,  carbon 
dioxide,  and  water  droplets  vary  with  altitude,  and  the  way  the  mean  - 
droplet  radius  varies  with  the  visibility.  As  better  values  for  these 
parameters  become  available,  the  model  can  be  easily  modified  to  ac¬ 
commodate  them. 

Tne  third  class  of  error  arises  from  basic  limitations  in  the 
modeling  process  itself.  Consider,  as  an  example,  the  transmission  of 
direct  solar  radiation.  If  we  use  best  estimates  of  the  manner  in 
which  the  concentrations  of  scattering  and  absorbing  particles  vary 
with  altitude,  we  can  predict  a  vertical  transmission  which  will  be  in 
error  by  some  amount  because  of  errors  we  described  as  class  two.  If 
now  the  sun  is  found  to  be  in  some  direction  other  than  vertical,  this 
error  will  be  amplified  because  of  the  greater  atmospheric  path  the 
radiation  traverses.  For  example,  if  our  vertical  transmission  error 
is  +10  percent,  a  solar  elevation  angle  of  30  degrees  results  in  an 
error  of  +21  percent  and  a  solar  elevation  angle  of  10  degrees  results 
in  an  error  of  +74  percent. 

Another  problem  with  a  model  of  this  nature  is  that  predictions 
of  properties  of  the  upper  atmosphere  are  being  made  based  on  ground 
level  measurements.  This,  of  necessity,  limits  the  cases  in  which 
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the  model  can  be  used.  For  example,  if  reflection  occurs  from  temper¬ 
ature  inversion  boundaries,  the  model  cannot  be  applied  to  cases  in 
which  an  inversion  layer  exists. 

Plans  for  the  next  six  months  include  several  phases.  First, 
some  limited  validation  measurements  will  be  made  in  the  visual  portion 
of  the  spectrum,  including  luminance  measurements  of  the  sun  and  the 
clear  daytime  sky.  These  measurements  will  be  converted  to  radiometric 
units  .aid  compared  to  predictions  made  with  the  model.  The  effect  of 
multiple  scattering  will  be  addressed  at  this  point. 

Second,  luminance  measurements  will  be  made  of  solar-irradiated 
armor  plate.  These  measurements  will  be  compared  with  predictions 
provided  by  the  model  after  using  available  bidirectional  reflectivity 
data. 

Third,  the  error  analysis  will  be  continued  in  an  attempt  to 
define  more  accurately  the  conditions  under  which  the  model  can  be 
utilized. 

Finally,  the  effect  of  the  ground  as  a  radiation  source  will  be 
considered  to  determine  those  cases  in  which  its  effect  may  be  ignored 
and  those  in  which  its  inclusion  is  mandatory. 
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Figure  1.  Format  Used  in  Predicting  Sky  Radiance 
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Figure  5.  6  as  a  Function  of  e.  and  ^  for  e-  3  60  Degrees. 


Function  of  e,  and  +  for  e  8  80  Degrees. 


Figure  7.  The  Scatte^ng  Area  Ratio  as  a  Function  of  2nR/x. 
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Figure  11.  Normalized  Directional  Mie  Scattering  Coefficient. 


Figure  12.  Altitude  Variation  of  the  Visual  Spectrum  Volume  Scattering  Coefficient 

after  WaldramH. 
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Mie  Transmission  Geometry  for  V  >  G(a). 
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Figure  15.  Carbon  Dioxide  Optical  Slant  Range  as  a  Function  of  9 
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Figure  16.  Carbon  Dioxide  Optical  Slant  Range  as  a  Function  of  e, 
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Carbon  Dioxide  Optical  Slant  Range  as  a  Function  of  $  , 


f'nncentration  as  a  Function  of  Altitude. 
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Figure  23.  Mie  Transmission  Geometry  for  V  <  G(x)  and  r  <  1/sin 


Transmission  Geometry  for  V  <  G(x)  and  r  >  1/sin 
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TABLE  I 


AX 

(Microns) 

HSCA) 

(10_3w/cm2) 

'em* 

Kr(X) 

o(X) 

YU) 

gA 

ST,R 

0.40-0.45 

8.40 

2.69xlOi0 

3.81 

0 

0.161 

0.0054 

0.0343 

0.45-0.50 

8.97 

1.70 

3.63 

0 

0.136 

0.0615 

0.0217 

0.50-0.55 

9.02 

1.13 

3.44 

0 

0.111 

0.3437 

0.0144 

0.55-0.60 

8.72 

7. 80xl09 

3.23 

0 

0.097 

0.4164 

0.0099 

0.60-0.65 

8.17 

5.55 

3.00 

0 

0.085 

0.1582 

0.0071 

0.65-0.70 

7.53 

4.06 

2.75 

0 

0.080 

0.0149 

0.0052 

0.70-0.75 

6.85 

3.05 

2.50 

0 

0.082 

0.75-0.80 

6.19 

2.33 

2.27 

0 

0.090 

0.80-0-85 

5.56 

1.81 

2.06 

0 

0.133 

0.85-0.90 

4.98 

1.43 

1.90 

0 

0.266 

0.90-0.95 

4.43 

1.14 

1.81 

0 

0.216 

0.95-1.00 

3.96 

9.24X108 

1.72 

0 

0.096 

1.00-1.05 

3.53 

7.59 

1.63 

0 

0.104 

1.05-1.10 

3.16 

6.22 

1.54 

0 

0.214 

1.10-1.15 

2.82 

5.21 

1.46 

0 

0.253 

1.15-1.20 

2.53 

4.36 

1.38 

0 

0.195 

1.20-1.25 

2.27 

3.70 

1.30 

0 

0.214 

1.25-1.30 

2.04 

3.15 

1.23 

0 

0.505 

1.30-1.35 

l.»5 

2.70 

1.16 

0.00214 

0.753 

1.35-1.40 

1.65 

2.33 

1.10 

0.00846 

0  829 

1.40-1.45 

1.50 

2.02 

1.04 

0.00779 

0.957 

1.45-1.50 

1.33 

1.76 

0.99 

0.00282 

1.106 

1.50-1.55 

1.24 

1.54 

0.94 

0.00282 

1.120 

1.55-1.60 

1.11 

1.35 

0.90 

0.00779 

1.044 

1.60-1.65 

1.03 

1.19 

0.87 

0.00846 

1.003 

1.65-1.70 

0.92 

1.05 

0.84 

0.00214 

0.904 
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